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We review the approach of generalized permutator to produce a class of integrable quantum 
Hamiltonians, as well as the technique of Sutherland species (SS) to map a subclass of it into 
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latter we describe how the complete spectrum can be gained by means of SS technique. Finally we 
speculate about possible applications to spin 5* models. 
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I. INTRODUCTION 



In recent years the discovery of materials that exhibit, 
at least in some energy regimes, a strong one-dimensional 
character has renewed the interest in models of interact- 
ing electrons in low dimensional lattices. 
Conducting polymers such as polyacetylene (CH)^ was 
probably the first noticeable example of such type of 
materials [1], where the interplay between 7r-electron 
correlation and chain dimerization was intensively stud- 
ied [2]. 

Successively, a very rich temperature-pressure phase di- 
agram has been observed for Bechgaard salts [3], which 
are linear-chain organic compounds such as (TMTSF)2X 
and (TMTTF)2X, where CIO4 or X = Br. In consider- 
able part of their phase diagram, namely at sufficiently 
high temperatures and pressure, the inter-chain hopping 
is irrelevant and such materials actually behave as ID 
systems. 

More recently, huge interest has been devoted to car- 
bon nanotubes [4] , which can be regarded to as graphite 
cylinders of a typical radius of some nm. Due to the 
discrete nature of the radial motion, only one band is ac- 
tually involved in determining the electronic properties 
(at least up to leV energy scales); for this reasons, as 
well as for the fact that they are ballistic, such materials 
are nowadays considered as the most promising material 
to probe electronic correlations in ID. 
Further significant example are semiconductor-based 
quantum wires synthesized through cleaved edge over- 
growth [5] , as well as Edge states in Fractional Quantum 
Hall Effect systems. 

In modelling any of the above materials, it is crucial to 
take into account the electron-electron interaction, since 
the latter heavily determines their physical behavior. It 
is indeed well known that in low dimension correlation 
effects become important, even for moderate interaction 



strength; this boils down to the general fact that in ID 
a single particle picture is not valid. The conventional 
techniques that are used in 3D to treat many body sys- 
tems are therefore either unreliable (such as mean field) 
or inapplicable (such as Fermi Liquid theory), when in- 
vestigating ID materials. For this reason, alternative 
approaches have been formulated and developed, both 
analytical (e.g. Luttinger Liquid Theory [6]) and numer- 
ical(e.g.DMRG [7]). 

In this general context, the availability of exact results 
for ID electronic models (at least for some values of the 
model parameters) plays a crucial role, not only as a 
conceptual outcome, but also as a comparison test for 
the above 'non-traditional' techniques. 

Within the exact-result approaches to one dimensional 
systems, the Coordinate Bethe Ansatz (CBA) is probably 
the most famous technique; it amounts to prove, when 
possible, that a given model Hamiltonian has eigenfunc- 
tions of the form proposed by Bethe [8] ; this method has 
been quite successfully applied to models of correlated 
electrons. For instance CBA has been used for solving 
the Hubbard model [9] and the t-J model [10], which are 
somehow the prototype models to account for correlation 
effects in electron systems. The former reads 

L 

^Hub = -^51 XI (SUS+I.'T + 4+1,<tS + 1,^) + 

L 

i=i 

where ^ and ^ are electron creation and annihila- 
tion operators on a ID lattice, j and a label position 
and spin of the electons on the chain, and fij^cr = Cj ^Cj ,^ 
is the number operator. In Eqs.(l) the first term is the 
tight-binding part and the second the on-site Coulomb re- 
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pulsion competing with it. The exact solution states that 
any small U > causes the ground state of the system 
to be insulating at half-filling, confirming the dramatic 
effect that correlations have in ID. 

The t-J model involves a spin-spin coupling term and its 
Hamiltonian reads 

L 

i ^ ^ 1 

where S = cj^^^y^c^-^, is the spin operator, and a —I 
(resp.t) if (7 =1 (resp.i). The exact solution of the t- 
J model for J = 2t was the first analytical proof that 
the Luttinger liquid hypothesis for 1-D systems was well 
grounded. Indeed the ground state of this model, for any 
filling value, is a liquid of singlet bound pairs of varying 
spatial separation and binding energy. The excitations 
were actually found to be of two types: spin- like and 
charge-like, and they can be proved to be gapless. 
The t-J model is also an example of constrained fermions, 
in the sense that the double occupied states are a priori 
excluded from the Hilbert space; this type of constraint 
physically emerges e.g. when investigating the low- 
energy singlet excitations of multi-band models; in the 
regime of strong correlation an effective 1-band Hamil- 
tonian can be obtained [11,12]. 

A quite detailed review of CBA-solved models can be 
found in [13]. 

Beyond the CBA, there exist another quite useful tech- 
nique within the framework of exact-result approaches. 

It is called Quantum Inverse Scattering Method (QISM) 
[14], and is in a sense more general than the CBA, in 
that it directly involves the notion of integrability of the 
Hamiltonian Ti. Here integrability actually consists in 
the possibility of finding a complete set of mutually com- 
muting observables 

[Jn,Jn']=0, (3) 

so that the eigenstates of H can be univoquely char- 
acterized by the quantum numbers related to such ob- 
servables. In the QISM, all the conserved quantities 
are encoded in an equation, known as the Yang-Baxter- 
Equation (YBE); each solution of the YBE determines 
a specific model, endowed with a set of conserved quan- 
tities. The standard way to determine the eigenstates 
within the QISM is again based on an Ansatz, called Al- 
gebraic Bethe Ansatz (ABA); differently from the CBA, 
the ABA explicitly lies on the underlying symmetries of 
the Hamiltonian contained in the related YBE. 

Within the above scenario, in this paper we first briefiy 
summarize in sec. II the main ideas of the QISM applied 
to fermionic models; then, in sec. Ill we review a recently 



developed method [15] to determine integrable electron 
models. This method, formulated within the QISM, 
amounts to finding solutions of YBE of polynomial form. 
We show in sec. IV that already the easiest case of first 
order polynomial yields a rich variety of integrable mod- 
els [16]. In section V we describe a technique (different 
from the usual ABA) to obtain the complete spectrum 
of the Hamiltonian for a subclass of the integrable cases. 
In particular, in sec. VI we review the integrable cases for 
extended Hubbard models (dimension of the local vector 
space dv = 4), whereas in sec. VII we explicitly derive 
the integrable cases for dy = 3, and look at their imple- 
mentation by means of constrained fermions. Finally, in 
sec. VIII we give some conclusions. 



II. INTEGRABILITY FOR ELECTRONIC 
MODELS 

A. Some introductory material 

In the second quantization formulation, a system of 
fermions on a ID lattice is described by c;reatioii and an- 
nihilation operators, which are governed by the algebra: 

{Ci,s, Cj,s'} = {Ci,s, c]s,} = Si J Ss,s' (4) 

where s = (a, a) is in general a multi-label, accounting 
for the orbital (=a) and for the spin (=(t). Here a can 
vary from 1 to A^orb and a can take 2 J -|- 1 values, where 
Noih is the number of orbitals and J the (half-integer) 
spin. 

In lattice model, one can associate to each site a local 
vector space Vj. As an example, in the case of s-orbital 
electrons we have a = s and a =t; ii ^ is made up of 
the 4 vectors 

Vj = Span[\^),=(^^\o),\D,=cl\o),; 

|o>.; |iT>, = c]^ct^|o),-) (5) 

where \o)j is the local vacuum. 

More generally, Vj has dimension dy and is spanned by 
vectors \a)j 

Vj = Span ^ \a}j = h'j^^ \o)j ,a = 1 . . . dy^ (6) 

Here the /ij"''s are products of creation operators cj^ 
with different s [17]. 

An important property of fermionic systems is that, 

due to the anti-commutation relations (4), the space 
Vj has an intrinsic graduation; this means that Vj can 
be decomposed into an odd and an even subspaces 
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^(1) 



V-^°^), where the odd (even) subspace 



Vj^"^ (Vj^^) is spanned by those vectors that are built 

with an odd (even) number of creation operators . For 
instance, in the case of the s-orbital electrons, the sub- 
space spanned by | "f) j and | J,)^ is odd, and the subspace 
spanned by \o)j and | is even. Similarly, the space of 
local linear operators acting on Vj (denoted by End(V^)) 
is also graded; odd and even vectors and operators are 
also said to be parity- homogeneous; in particular they 
arc respectively said to have parity p ~ I and p = 0, so 
that for any homogeneous 



e End{Vj) the rcla- 

^ , 15(05."^) +p(0f^) holds. Technically 

this is expressed by saying that End(y) is a graded local 
algebra. 



tion p(0(.") O^') 



Another important ingredient related to the vector 
space Vj, which plays an important role in the QISM 
to be presented below, is the state projectors. They are 
defined as 

Sji = \^)jM (7) 
The projectors fulfill very important properties: 



"7 "^Ja 



(8) 
(9) 



where [X, Y]± = XY - (-1)p(^)pW YX. 
In the case of s-orbitals they are explicitly given by the 
entries (a— th row and the 6— th column) of the following 
4x4 matrix 



4^(1 -n.^) 
(l-n^.^)(l-n.^) 

Si St 



S'tSi 



(10) 



t/ 



For constrained fermions the projector matrix is just 
the left upper 3x3 sub- matrix of (10). 

Notice that each of the entries of the above matrix is an 
homogeneous operator with parity p(fj^) — p{a)+p{(3) 
[18]. Projectors are useful in decomposing any opera- 
tor into its fundamental state-processes: any single-site 



operator O - 



(a) 



can indeed be written as 



of = (0W)^£:,-^ 



(11) 



where (O*-"^) is its representing matrix, defined through 
0(.'^)|/3), = (0W)^|a), (12) 



a functional equation for R. 

To each solution of the YBE one can associate an in- 
tegrable model; here below we shall briefly describe how 
to obtain integrable fermionic models (fermionic QISM, 
for more details see Ref. [19,20]). 

Prom the .R-matrix one can construct an operator- valued 
matrix Cj, related to each site j, as follows 



£|;(z.) = (-l)^(«)^'W^2?(«)^: 



■/3<5 * 
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-3^ 



(14) 



where the arc the projectors (7). Here a gives the 
row and /? the column of the entry of the Cj matrix (usu- 
ally called £-operator). Due to eqn.(8) the entries of two 
operators Cj and Ck of different sites fulfill 



B. Quantum Inverse Scattering Method for 
fermions: a short and schematic review 

The Quantum Inverse Scattering Method (QISM) is a 
powerful tool for studying quantum integrability because 
it provides models endowed with a set of mutually com- 
muting operators. Within the QISM a key role is played 
by the Yang-Baxter Equation (YBE) 



(I O R{u - v)) {R{u) O I) (I O R{v)) = 
= {R{v) O I) (I O R{u)) {R{u - t;) O I) 



(13) 



Here ^ is a rfy x C-number matrix, and I is the 

dy X dv identity matrix (dy being the dimension of the 
local vector space Vj); R depends on the complex num- 
ber u (called spectral parameter), and thus the YBE is 



Cj'^Ju),Ck]{u') 



ya,p,j,5 yu,u' yj^k 

(15) 



At the same time the property (9) of the projectors allows 
to show that the Yang-Baxter equation (13) is actually 
equivalent to 

R{u - v) {Cj{u) ®s Cj{v)) = (16) 

where the symbol 0^ is called ^rarferf tensor product, and 
it is defined as 

(A^ R-)({«}'{7}) _ _4{«}R{7}c_1^(K{«})-^K{/3}))K{7}) 



(17) 
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The presence of the above additional signs is strictly re- 
lated to the fcrmionic nature of the system, as observed 
in II A. 

Eq.(16) is a local relation, in that it only involves the 
j-th site of the ID lattice. However, using the property 



one can easily show that 



yjj^k (18) 



R{u-v){T{u)(3sT{v)) = (19) 



where 



T(m) ^ Cl{u) . . . C2{u)Ci{u) 



(20) 



Notice that T{u) is a operator-valued matrix (usually 
called monodromy matrix), whose entries are now oper- 
ators defined on the whole ID chain. 
Eq.(19) is the main relation of the QISM, because taking 
the trace of T(u) yields 



[trT{u),trTN{v)\ = 



Vm, V 



(21) 



which already implies the existence of an infinite num- 
ber of conservation laws; indeed developing the operator 
trT{u) in powers of ?i, eq.(21) states that all the coef- 
ficients (i.e. operators defined on the chain) mutually 
commute. 

Typically one searches for solutions of the YBE with the 
'boundary condition' that for a given point uq, R{u) is 
the identity matrix 



In this case one typically introduces 

Z{u) := {trT{uQ))-^trT{u) 

and defines 



Eq.(21) therefore implies 



n>l 







(22) 



(23) 



(24) 



(25) 



where Jn are the mutually conserved quantities. It can 
be shown that Jn is the sum of operators involving clus- 
ters of no more than n + 1 sites [21]. 
The first conserved quantity J\ is usually interpreted as 
the Hamiltonian, and the other ones as its symmetries. 
Ji has the form 



Ji = n 



(26) 



with periodic boimdary conditions 'Hll+i = T^Li- 
Each Ti-jj+i is a two-site Hamiltonian, and it is straight- 
forwardly connected with the .R-matrix solution of the 
YBE. Explicitly one can show that 



H 



(27) 



= (_l)p(7)(p(/3)+p(5)) duR"l{u) 



pip 5 



Writing down the 2-site Hamiltonian "Hj j+i in a ma- 
trix representation Hj j+i , the relation with the .R-matrix 
looks even simpler, namely 



(-^2 sites )-y^ 



(28) 



III. POLYNOMIAL R-MATRIX TECHNIQUE 

In the previous section we showed that to each solution 
R of the YBE (13) one can associate a fermionic Hamil- 
tonian 7i = Tij j+i endowed with a set of symmetries 
Jn?,. The derivative of R with respect to u basically 
represents the two-site Hamiltonian Hjj+i through the 
relations (40) (operatorial form) or (28) (matrix repre- 
sentation). 

The purpose is now to find solutions of the YBE, to 
figure out what kind of Hamiltonian is associated to the 
solution and to determine its physical features. Here be- 
low we describe a technique that we recently developed 
to find solutions of the YBE, and then we describe what 
kind of models we obtain. 

Notice that if R{u) is a solution of (13) then R' = 
f(u)R{u), with f{u) any scalar function satisfying 
/(O) = 1, is a solution as well, and it corresponds to 
an Hamiltonian H' = H + \1 where c = ^|„=o- 

Our approach consists in looking for solutions of (13) 

which are polynomials in the spectral parameter and sat- 
isfy (22) for uo = 0. This imphes: 



Riu) = I + W(F2 Sites + Cl) + l^^(') + . . . + ^^(f) 



(29) 



with p the degree of the polynomial (unknown a priori, 
and possibly infinite). Here {i = 1,. . . ,p) are un- 
known coefficients which are to be fixed through eqn.(13) 
and I the dy x dy identity matrix. 

Let us substitute the Ansatz (29) into eqn (13); due 
to the fact that YBE must be satisfied Vu, v, we have to 
equate all the coefficients of the same power u" u™. Thus 
we end up with algebraic equations, which can be grouped 
order by order, the order being the value I = m + n. 
The Zero-th and First Order equations are mere identi- 
ties which are always satisfied. For the Second Order we 
have 



^(2) = (^(1))2+ 



(30) 
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where = H + cl , and S, c are constants to be de- 
termined. 

The Third Order case still reduces to a single equation, 
which reads 

1^12'' ' ^23 } ^ ^23-^23 ~ '^-^23 ^12 ^23 + ^23^ = 

_rp(i)p(2)\ 6(1)6(2) 96(1)6(1)6(1)4.6(3) (0, 
~ 1 ^23 ' ^12 I -"-12 ^12 -^^12 -"'23 -"-12 ' ^12 ' V"' 

where { , } is the anticommutator. 

The successive orders in general consist of more equa- 
tions, up to the Highest Order 3p, in which case the 
equation is simply the (spectral parameter independent) 

YBE, 

f}{p) p(p) p(p) _ nip) f>(p) p(p) roo\ 
-'T-23 -'^12 -'^23 ~ -'^12 ^23 "-12 V>'^) 

Notice that, as the second order equation (30) is explic- 
itly solved once the Hamiltonian is given, the non trivial 
equations start from the third order; for a polynomial of 
degree p>\ one has 3p — 2 actual orders to consider. 
The idea is to proceed by successive attempts: once the 
Hamiltonian TL is given, one first tries to solve the equa- 
tions by means of a first order polynomial R -matrix; if 
this is not possible one passes to a second order poly- 
nomial and so on. If the Ansatz (29) is correct one can 
hope to find the .R-matrix in a finite number of steps, or 
guessing a recursive formula for the coefficient .R^'^ . The 
advantage of this method is that it deals with algebraic 
equations instead of functional equations (cfr.(13)). 

The simplest non-trivial case of a polynomial ii-matrix 
is of course obtained for p = 1. In this case R{u) = 
I + uM^\ and, due to (32), i?^-'^^ is nothing but a braid 
operator. Moreover, eq. (30) implies that the square of 
is proportional to the identity (since M"^^ = for a 
first order polynomial) 

{R^^^f oc I (33) 

Eqs.(32) and (33) imply that i?^^^ (and hence the 
Hamiltonian Ti.) is in fact an element of the symmetric 
group. The solutions of the symmetric group relations 
have been fully investigated for a two-dimensional local 
space (see for instance [22]); whereas the case of dimen- 
sion = 4 of the local vector space V was considered 
in [16], [20]. In the next section we extend the latter 
approach to generic dimension rfy. 

IV. GENERALIZED PERMUTATOR 
Ji-MATRICES 

It is well known that, among the solutions of the sym- 
metric group relations (32)- (33), one always finds the per- 
mutation operator P, defined by P{ea (X" e^) = e/3 (8) ea, 
where the e^'s for a = 1,. . . ,dv form an orthonormal 



basis for the local vector space V. In fact, this is true 
whatever the dimension of the representation is, and also 
holds for any graded permutation operator Pg, where 
Pg{ea ® ep) = (-)'("''(^)e^ ® e„, with e„ = 0, 1 grad- 
ing of the vector e^. Therefore, both P and Pg give rise 
to intcgrable models. One may wonder whether there are 
other solutions of symmetric group relations generalizing 
the structure of P, Pg. Here we propose the general- 
ized permutator ones, given by operators H which either 
transform one product of basis vectors into the reversed 
product, or leave it unchanged, according to a rule ex- 
plained below. In the matrix form this means that there 
is precisely one non-zero entry in each column and row of 
n. Moreover, due to (33), the non- vanishing entries of H 
must be equal to +1 or —1, up to an overall multiplicative 
constant. Explicitly, Va, f3 = 1, . . . rfy , 

n(ea 6^) = 6'^^(e„ ® e/3) -}- 6'°^(e/3(g)e„) . (34) 

diagonal tcrmts off-diagonal terms 

Here 9^^^ = Ca/js^^, and 61°^ = (1 - ea/s)^^^;' where e^^ 
is a discrete (0 or 1) valued function satisfying ea/S = e/Sa 
and €aa = 1) which selects the diagonal/off-diagonal 
terms; moreover, s^^ = ±1 accounts for additional signs 
of the diagonal entries, while = ±1 stands for the 
signs of off-diagonal terms; we shall impose = s^^ in 
order to make H a symmetric matrix. 

It can be easily verified that the operators H have 
square equal to the identity. Moreover, let us denote 
by S the set of index values, i.e. S = {1, . . . , dv}; the 
function ea/B on S x S characterizes the structure of any 
n. Let us denote by A"^ the subset of 5 x 5 of the couples 
{a, (3) for which ea/3 = 1: they determine the positions 
4 (a — 1) -I- /3 of the non- vanishing diagonal entries; since 
^afi = ^fia, whenever 6*^^ ^ we also have that 9^^ ^ 0; 
in other words, the subset A'^ can always be written in 
the form A'^ = \Ji SiX Si where the »Si's are disjoint sub- 
sets of S, with i = 1, . . . , Ns, and Ns < dy number of 
disjoint subsets in S . 

Inserting H into the relation (32) it can be easily found 
that it is a solution if and only if 6*^^ = Pi V (a, (3) G 
Si X Si, where pi = ±1 are signs. The values of the 
Pi's can be chosen independently, and the remaining off- 
diagonal non-vanishing entries are also free. In the fol- 
lowing we shall refer to these sohitions as the generalized 
permutators, and with abuse of notation we still denote 
them as H. 

Notice that this class of solutions exists for arbitrary dy, 
and in fact their actual number is an increasing func- 
tion of dv ■ As an example, in the next sections we shall 
investigate how these solutions provide 1440 integrable 
models in case of fermionic realizations with rfy = 4 (un- 
constrained fermions with two flavors) , and 56 models in 
case of two flavors constrained fermions {dy =3). At this 
stage, however, each of the integrable models proposed 
has no free parameter; we shall discuss in the next para- 
graphs how this limit can be partly removed by adding 
to the generalized permutators conserved quantities with 
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free parameters. Of course, the actual number of the lat- 
ter will depend again on the dimension of the local vector 
space. 

V. SUTHERLAND SPECIES AND SPECTRUM 

The meaning of a generalized permutator is easily un- 
derstood in terms of so called Sutherland species (SS). 
Starting from a given local vector space V, we may think 
to group its dv basis vector Ca (which in the sequel we 
shall denote as physical species PS) into Ng < dy differ- 
ent species Si which we call Sutherland species. Each of 
these species is left unchanged by the action of the gen- 
eralized permutator, the latter interchanging only basis 
vector belonging to different SS's. A generalized per- 
mutator would then have the structure of an ordinary 
permutator if represented on a local vector space of di- 
mension Ns- Notice that, with respect to this interpre- 
tative scheme, apparently the solutions presented in the 
previous section allow some extra freedom in the choice 
of signs of non-diagonal elements; on the contrary, here 
the off-diagonal elements connecting different PSs of the 
same SS Si to a different SS Sj must have the same sign, 
say Si. Such difference disappears through a redefinition 
of phase of some basis vectors belonging to Si ® Sj . 

The only signs in (34) which cannot be changed by a 
mere redefinition of basis are those of diagonal elements. 
Indeed, these signs can be used to classify the different 11 
solutions of GYBE for a given dimension dy- Following 
Sutherland's notation [23], each species can be classified 
as either 'fermionic'(F) or 'bosonic'(i?), according to the 
sign of diagonal elements. Hence all integrable models 
given in (34) are of type 

pNs-iQi , l = 0,...,Ns , 2<Ns<dv , (35) 

and characterize a model with Ns — I fermionic and I 
bosonic SS's. Different models that are recognized to be 
of the same SS-types also share the same structure of Co- 
ordinate or Algebraic Bethe Ansatz equations. The num- 
ber of different types of such structures associated with 
the generalized permutators can be easily evaluated, and 
is 

u{dv) = ^4 +'^dv-2 , (36) 

for given dy- We remark that this is not to be con- 
fused with the number of different integrable Hamiltoni- 
ans which in general is greater. 

The cases and BF^ have have been investigated 
in [23] within the framework of the Coordinate Bethe 
Ansatz, and ground state energy equations have been 
obtained. Moreover, the cases F^ and have been 
explicitly examined in [24] within the QISM; and some 
other cases {BF"^ and B^F"^) share the same algebraic 
structure as known models {t — J [25,26] and EKS [27] 



respectively). 

In particular, for the BF type of models, the entire 
spectrum can be obtained, and coincide, up to constant 
terms, with that of spinless fermions on an open chain: 

E{{ni}) =-2j2 cos (t^) m , (37) 
1=1 V / 

where ni = 0, 1 are quantum numbers and L the length 

of the chain. 

Notice that in the latter case we have the spectrum for 
any dy This means that in particular the ground state 
energy is independent of dy- However, it must be realized 
that the actual degeneracy of each eigenvalue depends on 
the way the Sutherland's species are realized in terms of 
physical species for each specific model. This enters ex- 
plicitly the calculation of both the ground state energy 
and the partition function, determining different physi- 
cal features of models with the same FB structure, as we 
shall discuss for the case = 4 in a subsequent section. 
Within the spectrum, the ground state is particularly 
worth of interest, because these Hamiltonians are ex- 
pected to well describe materials that exhibit peculiar 
physics at low temperatures. For the F^ case, it can 
be shown [23] that the miminum of the energy, whose 
eigenvalue reads 

is reached at equal densities of all fermionic species. 

Also, again as far as the ground state is concerned, it has 
been shown [28] that Sutherland's theorem (originally 
formulated for permutators of physical species) can be 
extended to the generalized permutators, and it is thus 
possible to assert that, in the thermodynamic limit, the 
ground state energy (per site) eo of a B"F"^ problem is 
equal to that of a BF'^ problem, for which the Bethe 
Ansatz equation have been formulated in full generality 
in [23]. 



VI. APPLICATIONS TO EXTENDED HUBBARD 
MODELS 

Here we review how the generalized pcrmiitator ap- 
proach succeeds in producing a class of integrable mod- 
els when applied to Hubbard-like systems of interacting 
electrons. 



A. Hamiltonian and representation 

We consider here a quite general 1-band extended 
isotropic Hubbard model preserving the total spin and 
number N of electrons, which reads 
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0',fe>,<T 



(39) 



In (39) t represents the hopping energy of the electrons, 
while the subsequent terms describe their Coulomb in- 
teraction energy in a narrow band approximation: U 
parametrizes the on-site repulsion, V the neighboring site 
charge interaction, X the bond-charge interaction, W the 
exchange term, and Y the pair-hopping term. Moreover, 
additional many-body coupling terms have been included 
in agreement with [29]: X correlates hopping with on- 
site occupation number, and P and Q describe three- 
and four-electron interactions. Finally, fi is the chemical 
potential. 



The local vector space V at each lattice site has dy = 4, 
and in the following we shall identify the 4 physical states 
I t), I i), |0) and I It) with the canonical basis Ca of C^. 

The matrix representation of the two-site Hamiltonian 
Hjj+i can be obtained through the projectors (10) in 
the following way 



w,,+i = (-i)^ww)+^w)(4^),)^;f 



Explicitly 
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where 



,22 
"22 



n + v-w 



h\l = hl\ = fi + V 



"is ~ "SI 



u23 _ l32 
"'23 ~ "■■ 



32 



M/2 



"21 



= w 



hlt = Y 



— /132 — —t 



h\\ = hli = t-2X + X 



"14 



"41 



hlt = htl = ^f, + P+^+2V-W 



"34 — "■43 



u 

^+2 



(42) 



hii = 2^^ + 4:P + Q + U + 4V-2W 



12 



7)1^ — h 
"34 — "43 



1^1 

'34 



-7)21 
"43 



t-x 



B. Integrable cases 

We would like to recognize in the matrix (41) one of the generalized permutators given in section 3. According 
to the scheme discussed there, this would guarantee that the corresponding Hamiltonian is then integrable. We first 

observe that this correspondence is possible only if t — X = 0; indeed only in this case the non vanishing entries of 

(2) 

^EH ™8,y coincide with those of a generalized pcrmutator. When this condition is implemented, the representation of 
the two-site extended Hubbard Hamiltonian (up to an additive constant c) is itself a generalized permutator whenever 
some linear relations among the non-vanishing Hamiltonian entries and the non-vanishing elements of a generalized 
permutator arc satisfied. It turns out that actually there are 96 different possible choices of values of the physical 
parameters in Ti satisfying such relations. They can be cast into six groups as (43) shows. 
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(43) 



Here s,- 



±1, 



1, . . . , 5 are arbitrary signs; the 



first two groups consist of 32 different solutions eacli, 
while any of the other four groups is made of 8 differ- 
ent cases. Notice that requiring the conservation of spin 
has reduced the 1440 generalized pcrmutators solutions 
of YBE for = 4 down to the present 96 solutions de- 
scribed in (43). 

As expected, t = X is a common feature exhibited by all 
the solutions, implying that the number of doubly occu- 
pied sites is a conserved quantity for those Ti that are 
derivable from first-degree polynomial i?-matrices. This 
feature is important in that it means that in these cases 
H can be diagonalized within a sector with a given num- 
ber of up and down electrons and doubly occupied sites. 
In practice, the solvability of the model in one dimension 
is not affected by having values of U and [i in Ti. other 
than those reported in (43) (see also [27]), as well as by 
the presence of an external magnetic field {i.e. adding to 
Ti. a term proportional to n^j — nj,|). 

The 96 integrable cases given in (43) can be classified 
according to the scheme of SS's, which allow to distin- 
guish just 12 different algebraic structures. In particular, 
the first group is characterized by Ns = dv{= 4) SS's, 
hence all models are of type F'^~^B^\ whereas the second 
and the third group correspond to Ns = 3 SS's, and the 
remaining three groups all correspond to Ns = 2 SS's. To 
all these structures are associated different sets of Betlie 
Ansatz equations, which have been explicitly worked out 
in [30], even though not explicitly solved. In what fol- 
lows we illustrate a different approach to some of these 
algebraic structures, which allow an explicit evaluation 
of many physical features. 



C. Ground state and thermodynamics 

Recognizing that a model identifies (up to some com- 
muting terms) a set of Sutherland Species greatly sim- 
plifies the calculation of the spectrum. The crucial point 
which allows that is the use of open boundary conditions, 
instead of the customary periodic ones; although in the 
thermodynamic limit the bulk properties are not affected 
by either choice, the calculations are more straightfor- 
ward for the former. Indeed in an open one-dimensional 
chain the set of eigenvalues of a generalized permutator 
is equal to that of an ordinary permutator between Ns 
objects, i.e. the effective dimensionality of the Hilbert 



space is reduced (reduction theorem). As a consequence 
of that, the degeneracy of the eigenvalues can also be 
computed, simply counting the ways one can realize a 
given configuration of Sutherland Species. 
The reduction theorem is proved when realizing that, ac- 
cording to what observed above, the relative order of any 
sequence of states belonging to the same species is pre- 
served [31]; the Hamiltonian can therefore be diagonal- 
ized within each subspace of given set of sequences. For 
instance, for Ns = 2 any model of FB type has a spec- 
trum given by (38) (plus conserved quantities), and the 
ground state energy per site is simply obtained by fill- 
ing all the lower energy levels up to the Fermi level; it 
turns out to be just a function of the total number of F 
particles in the ground state, np, and it reads 



eoinp) = 2nF sin{imp) + {U — f7)nt \+{iJ, — 



(44) 



here C7, and /x are the value of U and ji as fixed from ta- 
ble (43), n = ^ - fij, and n^^i = J2j "•i.T'^i.i itself 
a function of np, the explicit form of which depends on 
the specific FB model we are looking at. Hence, the true 
ground state energy of the model is obtained by mini- 
mizing (44) with respect to n^, and will depend on the 
physical parameter n and U. This is shown in fig (l)-(2). 



u 



10 
8- 
6 
4 

2- 


-2 
-4- 







*^ 






metal 


sula 




metal 


m-a 


c 




ni-b 


















II 





• (DO 



1 2 

n (filling) 



8 



FIG. 1. Ground state phase diagram of models with 
a single doubly degenerate SS of type F; for all models 
t — X — 1. Moreover the continuous lines describe the 
case X ^ {1 - a); Y = -a; P = -1; Q = 2 [28] ; the 
dashed line is the EKS model [27] {X = W = Y = V=1, 
P = Q = 0), and the dotted line corresponds to the AAS 
model {X = W = Y = V = P = Q = [31] ). All models 
exhibit an insulator-superconductor transition at n = 1, for 
different Uc- 

There we plot the ground state phase diagram for dif- 
ferent FB models, in which t = X = 1, U and n are 
varied, and the other interaction parameters are chosen 
in different specified ways, so as to have, within the same 
F species, two PS for (1), and three PS for (2). 
The empty circles represent empty states, the barred cir- 
cles represent singly occupied states, the full circles are 
doubly occupied states; regions in which there are no 
singly occupied states are characterized by tif = 0. 




n (filling) 

FIG. 2. Ground state phase diagram of the model 
X = X = l; Y = a; V = W = P = Q = 0. A filling 
controlled metal-insulator transition takes place for ?7 > 2 
between two finite regions, Ill-a and Ill-b. 

The dramatic differences in the phase diagram of mod- 
els represented in the figures are simply due to the fact 
that in the first case the F species contains two physi- 
cal species (| t), and | |)), implying n^^i = ^{np — n), 
and the B species contains the other two (| "fi), and |0)); 
whereas in the second case species F consists of one more 
physical species (| t|), hence n-f^^ = n — np), and B 
is just the empty state. Such difference also reflects on 
the degeneracy of each level of the spectrum, and ul- 
timately manifests in a different behavior of thermody- 
namical quantities. 

For instance, let us derive explicitly the partition func- 
tion for one of the models in fig. 1, corresponding to 
the Hubbard model with fixed bond-charge interaction 
{X = t) (also known as AAS model [31]). Implementing 



the relation between np and ?7||, the spectrum (37) in 
presence of Coulomb repulsion, chemical potential, and 
external magnetic field becomes 



E = E{{ni''^};Nr,Nu)= (45) 

k 



where = — 2tcosfc. 

The degeneracy g corresponds to the different ways one 
can realize a configuration of Sutherland species, with the 

constraint that the total numbers iV|| and appearing 
in (45) remain unchanged; a simple calculation yields 



g{E{{NF}-N^Y,N^)) = 



N 



IT 



(46) 



The rearrangement of the Fock space deriving from 

the idc;ntification of the Sutherland Species allows a 
straightforward calculation of the (gran-canonical) parti- 
tion function 



E E E 

L 



/(£;)e-'^^(^"n'^iT;^T) = (47) 

= (1 + e''^''-*))^ n + g[-/3(e*,-/^*(;*,/3,i7,fc))]^ 



k=l 



In the second line of (47) we have defined 



ti*{li,l3,U,h) 



1 



2 cosh (3h 



■exp2/3(/x- C//2) 



(48) 



/3 = l/lksT) being the inverse temperature. Notice also 
that the product over k resulting in (47) is in form sim- 
ilar to the partition function of a tight binding model of 
spinless fermions, where fi* plays the role of an effective 
chemical potential renormalized by the interaction U, the 
magnetic field h and the temperature itself. 

By means of the partition function (47), one can cal- 
culate the thermodynamic observables from the gran po- 
tential (per site) u = — liniL^oo ksT \n Z . This is done 
in [33]. In the following, we review the results concerning 
the specific heat. 

In fig. 3 we have plotted the specific heat (per site) Cy 
as a function of the temperature. Here the chemical po- 
tential n is eliminated in favor of the filling through the 
relation p = (n =)duj/dii, as usual. In particular, in 
the top figure we have examined the case of half filling 
(i.e. p = 1) and zero magnetic field {h = 0), for different 
values of the on-site Coulomb repulsion U. One can ob- 
serve that, across the value U/t = 4, the low-temperature 
behavior of Cy changes from linear to exponential; ex- 
plicitly, for U < At we have 



Cv ~ 7T 



with 



7 ^ 



6iV(i-(f/W) 



(49) 



whereas for U > At 



9 



{U - uf 



(i^-4t) 



(50) 



This is a finite-temperature effect of a metal-insulator 
transition, in accordance with the result obtained in [31], 
where a charge gap Ac = f/ — 4t is shown to open in the 
ground state for U > At. We recall that for X = (i.e. 
for the ordinary ID Hubbard model) no metal-insulator 
transition occurs; the bond charge term thus seems to 
give rise to a finite critical value Uc, increasing from to 
4:t as the coupling X is varied from to t. 
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FIG. 3. Tfio specific lieat as a function of temperature. 
Top: p = 1, h = 0: the metal-insulator transition is revealed 
tfirougli the change in low-temperature behavior from linear 
to exponential across the critical value U = At. 
Bottom: Cv for different filling values: p = 0.25 (dashed), 
p = 0.50 (dot-dashed); p = 0.75 (dotted) and p = 1 (solid); a 
low-temperature sharp peak emerges for non-vanishing mag- 
netic field. 

In the bottom fig. 3, Cy is plotted for different filling 
values, fixed ratio U/t = 2 and magnetic field h/t = 0.01. 
A sharp low-temperature peak, located at fcsT ^ h, is 
observed to emerge as soon as the magnetic field is turned 
on. 

One can show that 



lim hm Cy/T ^ lim hm Cy /T 



(51) 



differently from the ordinary Hubbard model, where the 
two limits are interchangeable [32] . At half-filling and for 
\U + \2h\ \ < 4t, for instance, one has Cy 7T with 



7 = 



6nt^(l-i{U + 2\h\)/4t)^) 



(52) 



Comparing eq.(52) to eq.(49), one can reahze that (51) 
holds. This can be interpreted as a signal of non Fermi 



liquid behavior. Similarly, the exponential behavior, oc- 
curring when the gap is open, is different; namely, for 
\U+\2h\ \ > 4t 



Cv ^ k] 



{U + 2\h\-Atf 



(C/ + 2|fe|-4t) 



{kBT/t)-^/'^e 5^ 



(53) 



to be compared to eq.(50). 

In fig.4 the specific heat of the Hubbard model with 

bond charge is plotted ior X = 1 is aside the case X — 
(i.e. the Hubbard model) for strong coupling, namely 
U = 8t. 




FIG. 4. The specific heat as a function of temperature for 
model the bond-charge model in the strong coupling regime 
{U = 8t), at half filling and zero magnetic field. The dotted 
line is the case X = - i.e. the ordinary Hubbaxd model- ob- 
tained from [16], and the solid line the case X = t, obtained 
from our exact calculations. Continuity arguments suggest 
that the specific heat for arbitrary < X < t lies between 
these two curves. The low-temperature peak originating from 
spin excitations is depleted by the bond-chaxge interaction. 

Notice that the ordinary Hubbard model has a low- 
temperature peak, whose origin is due to spin de- 
grees of freedom; the latter being not gapped, the low- 
temperature behavior of Cy is linear in spite of the fact 
that a charge gap is present at any U > [34]. In 
contrast, for X = t the spectrum does not carry any 
spin quantum number, due to the rich symmetry of the 
model; spins act therefore as a sort of dummy variables. 
Although the value X = t is a, particular one, it is rea- 
sonable to expect that, for continuity argument, the plot 
of Cv for intermediate values < X < 1 lies between the 
two curves. As a consequence, we can infer that the ef- 
fect of spin excitations is weakened by the presence of the 
bond-charge interaction, at least in the strong coupling 
regime. 
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VII. APPLICATIONS TO CONSTRAINED 
FERMIONS MODELS 

In this section wc would like to discuss a simpler ap- 
plication of the above scheme, which is the case of two 
flavors electrons (up and down spin, for instance) with 
the constraint that no two electrons can occupy the same 
site (constrained fermions). The local vector space in this 
case has dimension dy = 5, the physical species at each 
site being just | t), | i), and |0). 

The more general constrained fermions Hamiltonian with 

nearest neighbor interaction one may think of is an ex- 
tended t — J Hamiltonian [10], define as 

HeU = -t ^(4,aCj+l,<7 + h.C.) + jY^Sj- Sj+i 

j,<^ J 
+ Vj2njnj+i , (54) 



f)a 



(1 - nj,»)n 



■jat 



where Cja = (1 
and terms not conserving the total number of electrons 

and the total spin operator have been ignored. Interest- 
ingly, Heij reduces to the infinite U Hubbard model for 



J = V = Q, and to the standard t— J model for V = 
both of which have been solved in one dimension. 
Again, in order to answer the question whether there are 
HEtj Hamiltonians of generalized permutator type, we 
have to represent HEtj as a matrix, and this can be done 
by using the Hubbard projectors introduced in section 2 
(the constraint of no double occupancy identifies in fact 
the 3x3 upper left submatrix of (10)). Interestingly, it 
turns out that the representation of the 2 sites Hamilto- 
nian identified by (54) in this case has non- vanishing off- 
diagonal entries precisely where those of generalized per- 
mutators are. However, by solving the linear equations 
in the interaction parameters steaming from identifying 
such non- vanishing entries, it turns out that there are just 
8 different integrable Ti.EtJ, to be compared with the 56 
different generalized permutator obtained for dy = the 
reduction of solutions being due to the required conser- 
vation of spin and charge operators. 

The integrable cases can be classified according to the 
total number of corresponding SS into the following two 
classes: 



HEtA^l,S2) = -tJ2(^l,aCj+l,a + h.C.) + (si + Sa)* ^ fljUj+i 

j-o 3 

Heu{suS2) = -i ^(c]_^Cj-+i,a + h.C.) + 2sit ^(4" • Sj+i - -nj-fij+i) + (s2 + si)t^fijhj+i 

3,c 3 3 

(55) 



plus conserved quantities. Here the first group gives 
generalizations of the infinite U Hubbard model, the lat- 
ter being obtained for si = —S2\ whereas the second 
group generalizes the t — J Hamiltonian, which is in fact 
given by the choice si = l, S2 = — 1. 
Also, models in the second group all have algebraic struc- 
tures of type F^-i^i (these are four, in one to one cor- 
respondence with the four possible choices of signs ,si 
and S2); whereas models in the first group are of type 
p2-iQi particular, infinite U Hubbard model is of 
type FB (as expected from the known spinless fermion- 
like spectrum) , and t — J model turns out to be of type 

Interestingly, while integrability of models of the second 
group has already been discussed in [10], apparently 2 
{V = ±2t) out of the four models in the first group have 
not been studied elsewhere from this point of view. 
Among the models of first and second type the choices of 
Si = S2 are physically interesting, describing a model of 
interacting electrons in the limit of large on-site Coulomb 
repulsion and finite (= ±2t) neighboring sites Coulomb 
interaction. In the first group, for instance, the repul- 
sive case (si = S2 = 1) corresponds to a F^ model (the 
ground state and excited states of which being nicely dis- 
cussed in [23]), whereas the attractive case is a model. 
Both should be compared with the plane infinite U Hub- 
bard case (no Coulomb repulsion between neighboring 



sites), in order to understand whether neighboring sites 
Coulomb interaction can induce a quantum phase tran- 
sition. 



VIII. CONCLUSIONS 

In this paper we reviewed a recently developed method 
to determine integrable electron models, which amounts 
in finding solutions of the YBE of polynomial form. The 
method was applied to the simpler case of first degree 
polynomial, and was solved in full generality in this case 
(for any dimension of the on-site vector space); solutions 
are all and only generalized permutators. These were 
classified in terms of Sutherland species, which allow for 
a complete derivation of the spectrum in the case of mod- 
els of FB type, as well as for the calculation of the ground 
state energy in many other cases. We explicitly applied 
the scheme to the classification and solution of integrable 
electron models, both in the constrained and in the un- 
constrained case {dv = 3,4 respectively). 
To conclude, it must be stressed that in fact the scheme 
developed just requires to have an Hamiltonian which 
-when represented in matrix form on two neighboring 
sites- has the same structure of a generalized permuta- 
tor. Hence we expect the method to be capable of gener- 
ating integrable models also in physical context different 
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from correlated electrons. For instance, spin S Hamilto- 
nians have a local vector space of dimension dy = 25* + 1, 
implying that models discussed in terms of constrained 
fermions in the previous section could also have an in- 
terpretation as spin 1 Hamiltonian. Work is in progress 
along these lines. 
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